Exact results on the Kondo-lattice magnetic polaron 
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In this work we revise the theory of one electron in a ferromagnetically saturated local moment 
system interacting via a Kondo-like exchange interaction. The complete eigenstates for the finite 
lattice are derived. It is then shown, that parts of these states lose their norm in the limit of 
an infinite lattice. The correct (scattering) eigenstates are calculated in this limit. The time- 
dependent Schrodinger equation is solved for arbitrary initial conditions and the connection to the 
down-electron Green's function and the scattering states is worked out. A detailed analysis of the 
down-electron decay dynamics is given. 



I. INTRODUCTION 

The Kondo-lattice model (KLM) has found widespread 
application as a basic model for the description of itin- 
erant carriers interacting with local magnetic moments 
formed by inner shells of the constituting atoms. It 
was successfully applied in the theoretical description of 
the europium chalcogenides 1 " 3 -, aspects of the physics 
of manganites^— and the diluted magnetic (III,Mn)V 
semiconductors^. 

It has long been known that the KLM, although it 
forms a complicated many-body problem not solvable in 
general, has a non-trivial, exactly solvable limiting case 
- the (ferro-)magnetic polaroi d 10 ' 11 . This magnetic po- 
laron is set up by adding one electron into the other- 
wise empty conduction band and a ferromagnetically sat- 
urated background. 

Shastry and Mattisi^ gave a detailed discussion of the 
spectral properties of the polaron in terms of the re- 
tarded one-electron Green's function (GF). Berciu and 
Sawatzkyi 3 - have extended the GF solution to describe 
complex lattices and longer ranged exchange interactions. 
Knowledge of the GF allows for the direct calculation 
of the bare line shape of a photo emission experiment^ 
via the spectral density (SD). This strength of the GF 
approach comes at the expense of getting only indirect 
information about the underlying eigenstates of the one- 
electron quantum system. 

Therefore other approaches have solved Schrodinger's 
equation directly to get this state informatio n 1 ' 10 ! 12 ! 15 . 
Sigrist et. alJ£ gave a rigorous proof, that the ground 
state for a system with antiferromagnetic exchange cou- 
pling is of incomplete ferromagnetic order. We found 
these derivations to be incomplete in one way or another 
and there is, to our knowledge, no exhaustive derivation 
of the eigenstates for the finite system of N lattice sites 
in the literature until now. 

In the limit N — > oo the free dispersion of the electron 
becomes a continuous function of the wavevector. It was 
van Hove who pointed out in his seminal paper o 17 ' 18 that 
in case of a continuous (free) spectrum the interaction 
part of the Hamiltonian can lead to persistent effects 
not amenable to scattering theory. This self-energy ef- 
fects can be of dissipative nature (finite lifetime of quasi- 



particles) and/or cloud effects (formation of a new quasi- 
particle with infinite lifetime). Both effects are present 
in the problem of the magnetic polaron. 

The phenomenon of decaying states has attracted the 
interest of many physicists. Far from being comprehen- 
sive here we just want to mention the works, that have 
inspired us to the present investigation. 

Nakanishi 19 proposed an extended quantum theory, 
that contains complex eigenvalues of the Hamiltonian 
which can be associated with an unstable particle. To 
this aim he constructed the analytic continuation of the 
propagator for that particle and found poles in the con- 
tinuation. By introducing a complex distribution he was 
able to construct a new eigenstate of the Hamiltonian 
with complex eigenvalue (equal to the pole position) lead- 
ing to an exponential decay of the particle in time. De- 
viations from an exponential decay law at long times 
and the importance of the van Hove singularity on the 
lower bound of the spectrum was discussed by Hohler 20 
and Khalfin 2 !. The short-time deviations and the result- 
ing quantum Zeno effect where derived i n 22 ' 23 . Sudar- 
shan et. al. have extended the ideas of Nakanishi— and 
gave also a more formal derivation in a rigged Hilbert 
space formalis m 25 ! 26 . 

We will use the methods developed in the above men- 
tioned works as a pragmatic device to get a deeper un- 
derstanding of the decay dynamics of a down-electron in 
the problem of the magnetic polaron. 

The paper is organized as follows. 
In section [TTI we explain the model Hamiltonian and give 
the parameters used throughout this work. Section IIIII 
summarizes the results obtained with the GF approach. 
The spectral density of up/down-electrons is discussed 
in detail. The complete eigenstates for the finite lattice 
are derived in section [TV] Section [V] is devoted to the 
problems arising in the limit of an infinite lattice. New 
(scattering) eigenstates are constructed and the time- 
dependent Schrodinger equation is solved for arbitrary 
initial conditions. The connection of certain initial con- 
ditions with the electronic GF is given. In section I VII a 
detailed analysis of the dynamics of the quantum system 
for two different initial conditions is given. The connec- 
tion to the scattering states and the limitations of scat- 
tering theory are discussed. Finally we give a summary 
and draw conclusions in section l"VIII 
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II. MODEL 

Throughout the paper we are concerned with the s-f- 
like Hamiltonian: 



J=0.3eV 



H 



j- 



{z„Sfcfc irT + S°4_ a c ia ) (1) 



where a = ±, z± = ±1 and 5 4 = 5f ± iSf, describing 
free electrons hopping trough the lattice and undergo- 
ing a local (contact) interaction with immobile local mo- 
ments formed by the inner shells of the underlying atoms. 
The cf a {ci a ) denote the creation (annihilation) operator 
of an electron with spin a at lattice site R^. 
The above Hamiltonian constitutes a highly involved 
many-body problem not solvable for the general case. 
It has long been known, however, that there exists a 
non-trivial solution in the limiting case of one electron 
in the otherwise empty conduction band moving in a lat- 
tice of fully ferromagnetically ordered background spins 
(T = OK). 

The Hamiltonian ([I} does not distinguish energetically 
between different configurations of the localized spins 
in case of an empty band. One could enforce a fully 
aligned ground state by introducing an additional direct 
exchange term of Heisenberg type between the localized 
spins as was done in some worksi^. This additional term 
would lead to a true magnon dispersion of the localized 
moments (which is flat in our case) and to energy correc- 
tions of the eigenstates where a magnon is present. Since 
these energy corrections are small (typical magnon ener- 
gies are 2-3 orders of magnitude smaller than the elec- 
tronic energies) and the additional term would compli- 
cate the already involved calculations we omit this term 
and choose the fully aligned state out of the highly de- 
generate ground-state manifold "by hand". One has to 
keep in mind, however, that the dynamics of the localized 
spins are mediated only by the electron. 
If not stated otherwise we will use the following model- 
parameters for the concrete evaluation of our theory. We 
choose a magnetic moment S — ^ for the local spin sys- 
tem which reflects the magnetic moment of the prototypi- 
cal rare-earth compounds EuO and EuS. The underlying 
lattice will be the three-dimensional (3D) simple cubic 
(SC) lattice and the hopping Ty is chosen to give the 
tight-binding dispersion: 

W 

£k = — — (cos(fcx) + cos(ky) + cos(fcz)) (2) 
6 

with bandwidth W = 1. 



III. ONE-ELECTRON GREEN'S FUNCTION 

The limiting case of the "magnetic polaron" is usually 
discussed in terms of the retarded one-electron Green's 
function (GF). For the derivation of the latter one trans- 
forms ([T]) into reciprocal k-space and writes down the 




FIG. 1: (color online) Down electron spectral density along 
lines of high symmetry in k-space and quasi-particle density 
of states for two different coupling strengths J. Parameters: 
S = 3.5, W = 1.0 eV. 



equation of motion (EQM) of the GF. Under the above 
assumptions (ferromagnetic saturated local moment sys- 
tem, empty band) this EQM can be simplified in case of 
an up-electron (parallel to local moments) to give: 



1 



E + i0+ - e k 



JS 
2 



(3) 



An up-electron behaves essentially like a free electron 
with a shifted energy equal to the mean field value due to 
the ferromagnetic background. It has an infinite lifetime 
because it will not find a partner to flip its spin in the 
already saturated local moment system. 
The situation is quite different for a down-electron. In 
the EQM appears a higher spin-flip GF (SF-GF) and one 
has to go one step further in the EQM-hierarchy to get 
a closed set of equations. This results in: 



1 



E + i0+ - e k - £+(£) 



with the electronic self- energy^: 



. , JS f 
%l(E) = 1 



JG^(E- 



j&°)(E- 



(4) 



(5) 



where G^(E) = ^k^ '^) denotes the k-summed 
free electronic GF. The center of gravity of the down- 
spectral density is given by the first spectral moment: 



A I, 



(i) 



JS_ 



(6) 



which is the mean-field energy of the down-electron. De- 
spite this one expects significant changes in the down- 
spectral density compared to the mean-field result due 
to correlation effects. Especially possible spin-flips of 
the down-electron should result in a finite lifetime of 
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quasi-particles as is indicated by the complex- valued self- 
energy. Inspecting ([5]) reveals, that states with finite life- 
time can be expected exactly in the energy range of the 
up-spectrum because the self-energy has a finite imagi- 
nary part there. These states are commonly called the 
scattering states. In order to fulfill ([6]) for the center of 
gravity one then should find additional states (of infinite 
lifetime) above (J > 0) or below (J < 0) the scattering 
states as real roots of the denominator of (|4|) at least for 
sufficiently large JM- In the next section we will show, 
that there is exactly one such additional state called the 
bound state or polaron state. 

To illustrate our findings we have plotted the down- 
spectral density and the quasi-particle density of states 
(QDOS) in Fig. Q] for two different values of J. For large 
enough J (J — 0.3 eV, upper figure) the bound states 
are completely separated from the scattering states and 
form a polaron band. It can be show n 10 ! 11 , that in the 
limit of large J (\J\ WIS) this polaron band has its 
center of gravity at E? — ^(5+1) and the bandwidth is 
reduced by a factor of 2 s+i com P are d to W. 
When J is smaller (J = 0.1 eV, lower figure) parts of 
the polaron-band dip into the scattering states. There 
is still a visible polaron dispersion in the scattering re- 
gion, but now in form of a peak with finite linewidth 
implying a finite lifetime of the quasi-particle. We will 
show in a later section (|VI|) . that each quasi-particle peak 
in the scattering states is associated with a pole in the 
suitably analytically continued propagator leading to an 
exponential contribution in the decay dynamics. The po- 
laron dispersion is strongly disturbed at the positions of 
the van Hove singularities (horizontal dashed lines) of the 
scattering spectrum. The special role of the van Hove sin- 
gularities as limit /branch points for analytic continuation 
will also become clear in that section. 



IV. THE FINITE SYSTEM - EIGENSTATES 

In this section we derive the complete eigensystem for a 
finite lattice (N lattice sites) with periodic boundary con- 
ditions. Although parts of this derivation can be found in 
the literatur o 1 ' 12 , to our knowledge, the complete spec- 
trum was never calculated. Especially the appearance of 
pure up-electron states with one magnon emitted is not 
recognized. 

We use the following notation. The state of one electron 
with wavevector k, spin a and all local moments aligned 
(magnon vacuum) will be denoted by: 



-kCT 



|0;0) = |kcr;0>. 



(7) 



An up-electron with wavevector k plus a magnon of 
wavevector q is written as: 



1 



5-c+ |0;0) = |kt;q) 



(8) 



These states span the Hilbert subspace we are interested 
in. The Hamiltonian (JT|) commutes with the z-component 



of the total spin operator Sf ot = + \ S<r z o n ia)- 

Therefore we can classify the eigenstates by S* ot and by 
their (outer) wavevector k due to translational invari- 
ance. 

In the subspace of Sf ot — NS + \ the eigenstates are 
simply the up-electron states with magnon vacuum: 



ff|kt;0) = Uk 



-;0>. 



(9) 



The subspace S^ ot = NS — i is more interesting. It will 
be spanned by the states of one down electron in magnon 
vacuum and an up-electron plus one magnon emitted. By 
using the following Ansatz for the wave-function: 



\n)=A{a k \kl;0) 



q 



/3 k , q |k - q t; q) 



(10) 



we get a system of equations for the coefficients from 
Schrodinger's equation: 








This are N + 1 equations for the same number of co- 
efficients and we expect N + 1 eigenvalues per k-value 
(N(N+ 1) eigenvalues in total). The characteristic poly- 
nomial of (|llj) turns out to be of the simple form: 




with 



En 



E, 



E 

E 



~2~' 



JS 
~~2' 



(12) 

(13) 
(14) 



Not all e k - q „ will be different by symmetry arguments. 
Therefore we collect all equal e k -q„ in groups, where the 
numeric value of the group members is denoted by 

£k- Let us assume that 



with the convention that e 



(i) 



Nq such groups with respective degree of degeneracy g n 
exist. The q belonging to one group are denoted by qj 



with I G • • • g n — 1. With the definition: F n 
j_s 

2 



E- 



4^- we can recast (TT2|) to give: 




J Na 

3=1 3 



» 



(15) 



The first product contributes Nq differing eigenvalues: 



E 



{n) 



» 



JS_ 
~2~' 



(16) 
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The respective multiplicity of these eigenvalues will be 
g n — 1, since we divide by Fj (for j > 1) in the second 

factor of (fT5j) . Only E^p has the full multiplicity g\ . For 
the construction of the eigenvectors we notice first, that 
the eigenenergies (fTB| are equal to the pure up-electron 
energies ((9]) found in the Sf ot = NS + \ sector of the 
Hilbert space. It is therefore reasonable to assume the 
sought-after states to be up-states as well. We try «k = 
and /3k, q ^ only for q G {qj ^ } in our Ansatz (flO|) . 
From (|TT|) we then get the condition: 



E 

qefql"'} 



0. 



(17) 



This condition can be nontrivial fulfilled when the /3k, q 
are chosen to be suitable normalized powers of the prim- 
itive <?„th roots of unity. With this we find for the eigen- 
states: 



1 



/9r, 



g n -l 

1=0 



-ml I 



q| n) t;q, (B) > 



(18) 



where m <G 1 • • • g n — 1 is numerating the eigenstates of 
the subspace with energy E^ . 

One eigenstate is missing so far, since the degeneracy of 
is <7i. The states, which have only an up-electron 
(plus magnon) component are exhausted by (TT8"|) . By 

substituting £7^ into ([TT]) a simple calculation shows, 
that the last eigenstate is given by: 



|k|;0) 



■Ef 



1 



2iVS 
91 



(19) 



It has a finite down-electron component. 
Until now we have found 1 + (ft ~ = 1 + N ~ N G 
eigenvalues and corresponding eigenvectors from the first 
factor of (fT5j) . The remaining Nq eigenvalues come from 
the second factor: 



T N G 

F —X^El 

1 2N ^ Fj 

j=l 3 



(20) 




E 



= E 



The function on the rhs of (|2"(jj) has Nq — 1 zero- 
crossings between the up-state energies (|16p and one 
above (J > 0) or below (J < 0) the scattering states 
near E c s» ek + ^(S 1 + 1). The situation is depicted in 
Fig. [2j Rewriting (|20|) we get an implicit equation for the 
(polaron)-eigenenergies: 



E, 



with 



,p,(ra) 



JS 
~2 



*k(*) = 



1 E 



£k-q — £k 
2! - £k-q 



(21) 



(22) 




FIG. 2: (color online) Plot of eigenvalue equation (|20|l with 
positive (left figure) and negative (right figure) J. Parame- 
ters: |J| = 0.1, 5 = 3.5, e k = 0, 4 J) = -0-2, <4 2) = - al > 



J 3 ) 



0.1, 



„(4) 



0.2, 



The corresponding eigenstates are again obtained from 
(fTTj) to yield: 



|^' (n) ) = AA|k|;0) 



S q 1 



(23) 
|k-qt;q) 



with the obvious normalization factor: 



£k-q — £k 



P.O) 



Ck-q + — 



(24) 



THE INFINITE SYSTEM 



In the limit of an infinite lattice (N — > oo) the free 
electronic band energies ek become continuous functions 
of the wavevector. All summations will therefore be re- 
placed by integrations: 



^E-^/ dk 



(25) 



In this limit all eigenstates with down-electron part 
whose energies lie in the scattering region (ejj 1111 — 4jr < 
E < e£ lax — 4p) lose their norm as becomes clear im- 
mediately by inspecting (IT9l and (j24)l . Only the bound 
(polaron) state above (J > 0) or below (J < 0) the scat- 
tering region will survive provided that J is sufficiently 
large so that the state is energetically well separated from 
the scattering states. 

The question then becomes, what are the eigenstates in 
this limit? We will construct scattering states and show, 
that these states (+ polaron state) form a complete basis 
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in the continuum limit. Thereafter we solve the time- 
dependent Schrodinger equation which does not suffer 
from such difficulties for arbitrary initial conditions and 
show the connection to Green's function theory and the 
scattering states. 



A. scattering states 

In this section we ask for the result of a scattering 
process of an up-electron with a magnon, that is we want 
to solve the Lippmann-Schwinger equation^ : 

Ik-qtjq^ = |k-qt;q)+i?±i?i|k-qt;q) ± , (26) 

with the free propagator R± = (Eq — Hq ± i0 + )^ 1 . First 
we divide the Hamiltonian ((T|) into a free part Hq and an 
interaction part H x : 



Hn 



Hi 



£ ( e k - z a ^- | /z kiT . 



(27) 



J 



2VN 



kqcr 

JS. 

2~ 



With this division one can derive: 

iZ±fll|k-qt,q) = \Zi) + \Z 2 ) 
R±H X \Z X ) = A X \Z 2 ), 
R±H X \Z 2 ) = A 2 (\Z 1 ) + \Z 2 )), (28) 

where the Aj and \Zi) are given in the appendix (|A1I) . 
Defining the matrix: 



A =I^),UDU < 



(29) 



we can iterate (|26| and get: 



|k-qt> ± = |k-qt)+(j > ) (£a 



,t) + (;)u(gD.)u-.(jg 



|k-qt) 



1^2) 



1-A 2 (1 + A!) ^-^2 
= ^|k|;0)+^^ q ,|k-q't;q') (30) 



q' 



for the ingoing and outgoing scattering states. The diag- 
onal matrix D and the matrix of the eigenvectors U are 
given in the appendix lAl as equation (IA3|) . the coefficients 



aq and /3q iq / as equation (IA4|) . 

It remains to show, that the so constructed scattering 



states form a basis set in the limit of an infinite lattice. 
Applying the Hamiltonian results in: 



H\k 



£k-c 



JS 



|k-qt;q) ± (31) 



XX, q ,|k-q' t;q'>. 



q' 



The residue term 7£ q q , that is given in the appendix 
(|A5p goes to zero as 0(1/N). Therefore the scattering 
states (f30|) indeed become the correct eigenstates in the 
infinite lattice limit with densely lying eigenenergies in 
the expected region. 



B. time dependent Schrodinger equation 

Choosing the Ansatz: 

|*k(t))=Ok(*)|k|;0) + ^;^ k , q (*)|k-qt;q) (32) 

q 

for the wave-function one gets the following system of dif- 
ferential equations from the time-dependent Schrodinger 
equation (h = 1): 



id k (i) 



£k + a k (t) - Jsj A &,q(*)> 

q 

S_ 

2N~ 



£k-q 
J 

"2N 



JS\ 



£&,,(*)■ 



(33) 



This can be transformed into a system of algebraic equa- 
tions by the Laplace transform: 



e~ s 7(t)dt; Re(s) > 



(34) 



and one obtains the solution for the coefficients in the 
s-domain: 



<*k(s) 



«a k (0) 



US 



JS 

2 



JG<°)(is+^) 



fl 1 JG 

y 1 i_i G (o) (is+ ^> ) 

7 Ac,q(0) 

\'2.YS^ (is-e k _ q + 4£)/n(is) : 



(35) 



0k,q(s) 



*j8k,q(0) 



ia k (0)-(is-ek + j y )ok(») 



^ - £k-q + 4f 



with 



fel(2r) = «- Ck-^ffir+l + <S>k(z + ^-) 



(36) 



(37) 



For the back-transformation into the time-domain one 
has to solve the Bromwich integral: 

e st f(s)ds, (38) 
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FIG. 3: (color online) Integration contour in complex plane. 



where 7 is a real constant, that is larger than the real 
part of any singularity of f(s). In our case all singu- 
larities/branch cuts of ak (s) and /?k,q(s) he in a finite 
domain on the imaginary axis. Therefore 7 can take any 
value 7 > 0. The integration path is shown in Fig. [3] and 
denoted by G\. The contour can be closed at infinity 
since the integrand vanishes there (dashed black line). 
Performing the variable substitution s — > — iz and a con- 
tour deformation one can change the integration path to 
become C2 (red line). The singularities/ branch cuts now 
lie on the real axis. 

We will give the solutions for two special boundary con- 
ditions here. The choice: 



a) Ok(t = 0) = 1; /3 k , q (t = 0) = 0, Vq 



(39) 



results in: 



aik(t) = — ^ / (1: 

ZTTI 



C2 z- 

OO 



c .75 f-i , JGW(z+4f) 
£k 2 \ L T 1- J G (°)(z+4f ) 

— / dx e-«* (G u (x + iQ+) - G n (x - »0+)) 

/OO 
dxe~ txt S ki (x), (40) 
-00 

that is, ak(t) is the Fourier transform of the spectral 
density obtained from the down-electron Green's func- 
tion derived earlier Q. 
Choosing: 

b) /3k, q „(0) = l;ak(0) = /3k, qm (0) = 0,Vm ^ n (41) 



one gets: 



«k(i) = 



JS 1 



(42) 
of the 

scattering states (I3U1) for large times. We will test this 
conjecture in the next section, when we have build up the 
machinery to investigate different time domains of a^(t). 



It is to be expected, that this will become the 
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FIG. 4: (color online) Probability to find a down-electron 
P(t) — \ak(t)\ 2 as a function of time and the free electron 
dispersion ek for two different values of J. The time is given 



in units of ; 
eV. 



6.58212* 10 _lb s. Parameters: S = 3.5, W = 1.0 



VI. DYNAMICS OF DECAY 

The absolute square of the coefficient a-^(t) in (|32|) can 
be interpreted as the probability to find a down-electron 
in the system at time t. In Fig.@]we show P(t) — |ak(i)| 2 
of a system subjected to the boundary conditions (|39|) for 
two different couplings J = 0.1/0.3 eV. For large J (up- 
per panel in Fig. Q] shows the corresponding spectral den- 
sity) there is always a finite probability to find a down- 
electron in the system. P(t) shows characteristic oscilla- 
tions over time. These oscillation will be damped with 
increasing t and P(t) becomes static in the long-run limit. 
When J is smaller, parts of the polaron band dip into the 
scattering states (lower panel in Fig.[T]). For wavevectors 
k where £k lies below a certain threshold P(t) decreases 
exponentially and become zero in the long-run limit. The 
down-electron decays inevitably. However, it is clearly 
visible in the figure, that there are deviations from an 
exponential decay law in the short and the long time 
regime. We will now give a more detailed analysis of the 
different time domains. 



A. short time behavior 



7 

Ci 



For small times a Taylor expansion of the exponential 
function in (|4T))) gives: 



OCk(t) 



1 + ixt 



(xt) 2 .(xt) 1 



G 



A/, 



(2) 

ki t 2 



A/, 



(3) 

ki ,3 



The spectral moments , can be calculated exactly in 
our case and we give the first four in the appendix (fCl) . 
From this we get the short time expansion of P{t): 



G 



^ + 0(t 4 ). (43) 



|a k (t)| 2 = l-dt 2 + 0(i 4 ) 



with 



Ci = M, 



(2) 



(M^ 2 



J 2 5 



(44) 



(45) 



Since Ci > P(t) will always decrease (or stay con- 
stant) for small times as it should be for normalization 
reasons. The change of P(t) has a zero slope at t = 
since dP(t)/dt\ = 0. This leads to the famous quantum 
Zeno effect— - a down electron whose existence is moni- 
tored continuously by measurement— will never decay. 



B. intermediate and long time behavior 

For times t > we can omit the part of the integration 
contour C2 in (|4U)l that lies below the real axis because it 
does not contribute to the integral. The remaining inte- 
gration path is shown in Fig. [5] as C\. To study the long 
time behavior we pull the contour into the lower com- 
plex plane. For this aim we mention first, that the down 
electron GF has a branch cut in the energy region of the 
scattering states. To pull the contour below the real axis 
we have to find the analytic continuation of the GF over 
the cut from above to below the axis. This is equivalent 
to find the analytic continuation of the free GF: (z) 
as becomes clear by inspecting (j40|) . To be specific we 
choose here the free GF of the simple cubic lattice with 
the dispersion given in @. The analytic continuation 
of this function is discussed in detail in the appendix 
(|B|) . The important point is, that the van Hove singulari- 
ties are branch points and one reaches different Riemann 
sheets depending on the position relative to the singular- 
ity (left/right) where the continuation is performed. 
The situation is depicted in Fig. [5] The green dots mark 
the positions of the van Hove singularities. The deformed 
contour C2 (red) cannot overcome these points analyti- 
cally. The contributions to the integral (|40l) from the 
path left and right of the singularities will not cancel 
each other because they lie on different sheets of the Rie- 
mann surface. We exemplify the (approximate) analy- 
sis of the integral by giving a concrete expression for the 




FIG. 5: (color online) Deformation of the integration contour 
in the complex plane to study the intermediate and long time 
behavior. Green dots: van Hove singularities of the scatter- 
ing states. Yellow dot: separated polaron pole. Black dot: 
complex pole of decaying state. 



path around the lower band edge at x s 



1/11 



11 



-1 



— OO 
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(46) 



2tt 



dx e 



-% _4„tjn t ( 1 



20? 



(lcl'H-i + o{t-i) 



(47) 



The other van Hove singularities give rise to similar terms 
and we have summarized the explicit expressions for the 
coefficients C s in the appendix (|D|). From (l47l) we see, 
that the van Hove singularities contribute an oscillatory 
term with a frequency equal to the energetic position to 
ak(t). This term is damped in time with a decay rate 
given by a power law. The exponent of this power law 
depends only on the nature of the van Hove singularity, 
that is a square root singularity in our case— 
For large enough J there is a separated pole of the down- 
electron GF on the real axis shown as yellow dot in Fig. [5] 
If we denote the position of this pole by x po \ the contri- 
bution to a k (t) is given as the residue of the integrand 
in (j40|) at this point: 



a 



pol 



(t) 



1 - JSd x 



Cpol^ 



-ix vo \t 



(48) 



The characteristic oscillations of P(t) visible in the upper 
panel of Fig. [4] are superpositions of the oscillatory terms 
in (j47]l and (|48| . They are damped in time as (~ t~*) 
at least. For large times only the contribution from (|4"8")l 
will survive and we get: 



lim |« k (t)| 2 = |C pol | 



(49) 



For large coupling strength |J| ^> 1 the pole position 
becomes approximately x po i ss ^(S + 1). Using the high 
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energy expansion of the free GF in (l48l) we can derive an 
explicit result for P(t) in this limit: 



&5,M*>I * 25TT 



(50) 



The electron polarization can be calculated from P(t): 



(**) (*) = 



EJ/W*)| 2 -M*)| S 



|a k (i)| 2 , (51) 



where in the last step the normalization condition 
E q IAc, q M| 2 + M0| 2 = lis used. 
We come now to the case of small \J\, when parts of 
the polaron band lie inside the scattering states. For 
those k-values the initial down-electron will decay over 
time and vanish completely in the long time limit. One 
therefore expects an exponential decaying contribution 
to ctk(t). As soon as the polaron peak runs into the scat- 
tering states there appear poles in regions II, III or IV of 
the analytic continuation of the down-electron GF that 
account for such an exponential term. This is depicted 
in Fig. [5] by the black dot in region IV. These poles con- 
tribute a term to ak(t) that equals (|48|) but now with a 
complex pole energy a; po i = x r — £7. 
There is one pole in each sheet of the different analytic 
continuations. In Fig. [6] (lower panel) the pole lines in 
the complex plane (parametrized by 6k) are shown for 
J = 0.1 eV. Whenever a pole has a real part x r , that lies 
inside the energy region where the analytic continuation 
of the corresponding sheet has been performed (II, III or 
IV) it will add a term (|48|) to ak(t). The upper panel 
of Fig. [5] shows the free dispersion ek for a path along 
the high symmetry lines in the first Brillouin zone. The 
6k for which the corresponding pole is inside its region of 
visibility can be read off from the horizontal colored lines. 
Whenever this is the pronounced (quasi particle) 

peak is visible in the down-electron spectral density, that 
lies approximately at the position of the real part of the 
pole x r . We demonstrate this in Fig. [7] for two different 
values of ek- For ek = —0.5 eV there is only a pole in 
region II (endpoint of the red pole-line) . One finds a dis- 
tinctive quasi-particle peak in the spectral density (red 
line) at approximately the real part of the pole position 
(vertical dotted line). An interesting constellation arises 
for 6k = —0.23 eV. From Fig. [S] one can read off that two 
poles are now in their respective region (III, IV). Con- 
sequently two quasi-particle peaks can be found in the 
SD (black line) . The negative imaginary part of the pole 
position can be interpreted as (half of) the inverse life- 
time: r ~ ^. This is reflected directly by the width of 
the quasi-particle peaks, the smaller r the broader the 
corresponding peak in the SD. 

We had seen, that there are deviations from an exponen- 
tial decay law at small times. The same is true for the 
long time limit. When t > ^ the exponentially decaying 
term 1(15]) is damped away and only the contribution from 
the van Hove singularities (H7|) remain, which leads to a 
power law decay rate (|ak(£)| 2 ~ i~ 3 )- It was first shown 
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FIG. 6: (color online) Lower figure: Pole lines (parametrized 
by 6k) in the different Riemann sheets. The colors of the lines 
denote the sheet where the poles can be found: red: II, green: 
III, blue: IV. When a pole lies in its window of visibility 
(dashed area with the same color) it gives an exponentially 
declining contribution to Qk(i). The yellow line shows the 
positions of the polaron peak, when it lies outside the scat- 
tering region. Upper figure: Free electron dispersion ek and 
the windows of visibility for the different poles. Parameters: 
S = 3.5, J = 0.1 eV and W = 1.0 eV. 




FIG. 7: (color online) Down-electron spectral density for two 
different values of tk- Vertical black, dashed lines: position 
of the van Hove singularities. Vertical dotted lines: real part 
of the pole positions. Parameters: S = 3.5, J = 0.1 eV and 
W = 1.0 eV. 



by Khalfin 21 , that this deviation from exponential decay 
at large times is a general property of a system with a 
spectrum, that is bounded from below. 
At the end of this section we want to discuss the long 
time behavior of a system, where we have chosen the 
initial conditions (1411) . If we deform the integration con- 
tour C2 in (|4"2"1) as shown in Fig.[5]and chose J sufficiently 
small, so that h\(z) has no real root outside the scatter- 
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ing spectrum we get for large times: 



lim a k (i) = c 

t— foo 



-i(ck- 



2 ) h s 



(52) 



V27VS e k - q - e k - ^(25 + 1) - ^$ k (e k 



z0+) 



This is indeed identical to the coefficient found earlier 
for the scattering states (|30]) up to a trivial time depen- 
dent factor. For larger J there will be additional terms 
coming from the real root of h\{z) which are not con- 
tained in the result of scattering theory. The reason for 
this are the persistent self-energy effects (cloud effects) 
already mentioned by van Hoveil. 



VII. SUMMARY AND CONCLUSIONS 

In this work we have revised the theory of the magnetic 
polaron to get a deeper understanding of one of the few 
exactly solvable limiting cases of the KLM. 

First we have derived the complete eigenvalue spec- 
trum of the finite lattice and the corresponding eigen- 
states are constructed. 

In the limit of an infinite lattice parts of these eigen- 
states lose their norm and the eigenvalues of these states 
degenerate with the eigenvalues of the pure up-states. 
We then calculate the scattering states and show, that 
these states become the new eigenstates in this limit. 

By solving the time-dependent Schrodinger equation 
for arbitrary initial conditions we are able to give a de- 
tailed analysis of the decay dynamics of the system. This 
is done for a down-electron with wavevector k prepared 
at t = 0. For large exchange coupling J the probability 
P(t) to find a down-electron at time t > shows char- 
acteristic oscillations which are damped as ~ t~?. We 
find the reason of these oscillations to be interference of 
different oscillatory terms with frequencies equal to the 
energetic positions of the van Hove singularities and the 
separated pole (polaron peak). For large times P(t) be- 
comes static with a value P(t) < 1 and for |J| 3> 1 an 
explicit value can be derived. 

For small coupling J and e k below a certain threshold 
P(t) goes to zero over time - the down-electron decays. 
The main contribution to the decay dynamics stems from 
poles in the analytic continuation of the propagator giv- 
ing rise to an exponential decay law. The imaginary part 
of the poles determines the lifetime of the down-electron. 

We find deviations from an exponential decay law for 
small and large times. At t = the slope of P(t) is 
zero. This leads to the quantum Zeno effect. At large 
times the exponential term is damped away and the decay 
behavior is solely determined by the contributions from 
the van Hove singularities. We get a power law decay 
(P(t) ~ t~ 3 ) in this regime and there are also oscillations 
caused by interference effects. 

The connection between certain initial conditions for 
the solution of the time-dependent Schrodinger equation 



and the retarded Green's function approach as well as 
scattering theory is worked out. When a pure down- 
electron state is prepared at t = the time development 
of the down state coefficient of the wave function is given 
by the time-dependent down-electron spectral density, 
that can be obtained from the retarded Green's func- 
tion. Preparing an up-electron and a magnon at t = 
and taking the t — ¥ oo limit the result of scattering the- 
ory (for the down-electron coefficient) can be reproduced 
for small J. For large J however a new quasi-particle 
with infinite lifetime is created. Scattering theory is not 
able to describe this formation of a new quasi-particle. 

The limiting case of the magnetic polaron is believed 
to describe the essential physics of the magnetic semicon- 
ductors EuO and EuS. For a realistic description of the 
quasi particle density of states of these materials one has 
to replace the model band structure (|2|) used in this work 
by a material specific one, that can be obtained from a ab 
initio band structure calculation. Such a calculation was 
done by Nolting et. al<2& for EuO and Borstel et. al^for 
EuS. From the experimental side the natural way of mea- 
suring the unoccupied states would be an inverse photoc- 
mission experiment. An alternative method is the two- 
photon photoemission spectroscopy^. This experimental 
technique can be used to get time and spin resolved spec- 
tral information with time resolutions down to the fem- 
tosecond and even attosecond regime^—. These new 
experimental methods could open up a possible road to 
a direct measurement of the decay dynamics of an (down) 
electron in the above mentioned prototypical compounds. 



Appendix A: Definitions scattering section 
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FIG. 8: (color online) Free GF slightly above the real axis. 
Vertical dashed lines: position of the van Hove singularities 
that form the boundaries for a possible analytic continuation 
below the real axis. Parameters: W = 1.0 eV. 
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Appendix B: Analytic continuation of G^(z 



clearly visible van Hove singularities at the band edges 
(—3,3) and at (—1,1) are limiting points for a possible 
analytic continuation over the branch cut. Therefore one 
reaches different Riemann sheets, depending on the sec- 
tion (II, III or IV) where the analytic continuation is 
performed. 

To find these continuations we discuss firstly the analytic 
properties of the constituting functions. The square root 
function has branch points (of first order) at 0, — oo and 
a branch cut connecting them on the negative real axis. 
Whenever one crosses this cut, one has to change the sign 
of the square root function to get an analytic continua- 
tion. 

The hypergeometric function: 



TT.vplci =_.F, ( l^Alx 



(B3) 



has a branch point of second order at the points 1, oo and 
a branch cut on the real axis connecting them. Using 
formulas given in (35j one finds for the continuation from 
above to below the real axis: 



Hyp2(z) = (-z)~i . 2 F X 



112 1 



+ (-*)" 



r 



r(l)r ~V3'3'3'* 
F (1 2.4 1 

2 2-fl [ 77, n > o) - 



(B4) 



and from below to above the real axis: 



Hy P 3(z) = { -z) ^ ( 



(r(D) 



+ {-*)- 



r 



(r(I))' 



3' 3' 3' z 



2 2 4 1 



The simple cubic lattice Green's function can be ex- 
pressed in terms of the Gauss hypergeometric function 34 : 



G<°>(z)=/ii(*) 



with 



2-^1 



1 2 1 ( \ 



(Bl) 



With this analysis it is now easy to find the analytic 
continuation of (jBlj) from above to below the branch cut. 
To exemplify this, we discuss the continuation in region 
II. By crossing the real axis one also crosses the branch 

cut of the square root function \J~\ + ^. Defining new 

functions: 



2z 



*\i i -i-Ji-i 



1 / . o . ,„ . 2s /- 9 „_ L 1 



(B2) V2(z) 




m{z) = ±z 2 - (9 - 4z 



(B6) 



This function is analytic in the complete complex plane 
with the exception of a branch cut on the real axis in 
the region E B <E -3... 3 (bandwidth: W = 6). When 
the real axis is crossed in this region, the imaginary part 
of (IB1I) shows a discontinuity of A.(x) = ±2irp(x) where 
p(x) denotes the free density of states at point x. In Fig. [5] 
G(°\x + i0 + ) is plotted slightly above the real axis. The 



one then finds: 



G^(z)=p 2 (z) (HyplMF))) 2 



(B7) 



for the analytic continuation. 

Similar results can be obtained for the other regions. 
They are summarized in Fig. |H1 The two missing defini- 
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They can be calculated algebraically by the following 
rulei^: 



[[...[[c kl) H]_ > H]_... > HUc+} + ). (C2) 



FIG. 9: Analytic continuations of the free GF from above to We give the first four moments for the limiting case of 
below the real axis over the different branch cuts. Parameters: ^} ie magnetic polaron' 
W = 1.0 eV. 



tions are: 
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Appendix C: spectral moments 

The moments of the down-electron spectral density are 
defined by: 
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Appendix D: coefficients long time expansion 
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A short summary of how the spectral moments are defined 
and how they can be calculated is given in appendix (fC| . 
A sufficient condition is, that the center of gravity lies out- 
side the scattering states, that is \ JS\ > W/2. 
"Continuously" means in sufficient short time intervals so 
that Cit 2 < 1. 



